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ABSTRACT 

Y CVn is a carbon star surrounded by a detached dust shell that has been imaged 
by the Infrared Space Observatory at 90 fmx. With the Nancay Radio Telescope we 
have studied the gaseous counterpart in the 21-cm Hi emission line. New data have 
been acquired and allow to improve the signal to noise ratio on this line. The high 
spectral resolution line profiles obtained at the position of the star and at several off- 
set positions set strong constraints on the gas temperature and kinematics within the 
detached shell; the bulk of the material should be at ~ 100-200 K and in expansion 
at ~ 1-2 kms . In addition, the line profile at the central position shows a quasi- 
rectangular pedestal that traces an 8 kms -1 outflow of ~ 1.0 10~ 7 M Q yr _1 , stable 
for about 2 10 4 years, which corresponds to the central outflow already studied with 
CO rotational lines. 

We present a model in which the detached shell results from the slowing-down of 
the stellar wind by surrounding matter. The inner radius corresponds to the location 
where the stellar outflow is abruptly slowed down from ~ 8 kms -1 to 2 kms -1 (ter- 
mination shock). The outer radius corresponds to the location where external matter 
is compressed by the expanding shell (bow shock). In this model the mass loss rate of 
YCVn has been set constant, at the same level of 1.0 10 -7 M^yr -1 , for ~ 4.5 10 5 
years. The gas temperature varies from ~ 1800 K at the inner limit to 165 K at the 
interface between circumstellar matter and external matter. 

Our modelling shows that the presence of a detached shell around an AGB star may 
not mean that a drastic reduction of the mass loss rate has occurred in the past. The 
inner radius of such a shell might only be the effect of a termination shock rather than 
of an interruption of the mass loss process. 

Key words: stars: AGB and post-AGB - stars: carbon - (stars:) circumstellar matter 
- stars: individual: Y CVn - stars: mass-loss - radio lines: stars. 



1 INTRODUCTION 

The history of the mass loss experienced by stars on the 
Asymptotic Giant Branch (AGB) is a key issue for describ- 
ing the late stages of evolution of low and intermediate mass 
(1-6 Mq) stars, as well as an important ingredient for the 
characterization of the cosmic cycle of matter. Determining 
the mass loss rates of red giants is generally based on mod- 
elling radio molecular line profiles (Schoier 2007) or infrared 
continuum energy distributions (van Loon 2007). However 
these methods are limited to the central parts of circum- 
stellar shells, and, as the mass loss rates of AGB stars are 
variable, it has been difficult to establish a balance of the 
mass loss over the long periods of time that need to be con- 
sidered (10 4 -10 6 years). 

In principle the H I line at 21 cm should be a useful 
tracer of AGB circumstellar environments (Le Bertre et al. 



2005). Hydrogen dominates their composition (~ 70 % in 
mass) and should be in atomic form, at least in the external 
parts of these shells (r ^ 0.1 pc). However the 21 cm line is 
weak and generally contaminated by the much more intense 
galactic emission arising on the same lines of sight. Hi line 
observations of AGB stars have long been limited to Mira 
(Bowers & Knapp 1988). However, since 2001, using the up- 
graded Nancay Radiotelescope (NRT), we readdressed this 
issue and detected several objects in emission, in particu- 
lar RS Cnc (Gerard & Le Bertre 2003), EP Aqr & Y CVn 
(Le Bertre & Gerard 2004, Paper I), and X Her (Gardan et 
al. 2006, Paper II). Recently, we also presented the results 
of a survey of 22 sources, with 18 new detections (Gerard 
& Le Bertre 2006, Paper III). Using the Very Large Array 
(VLA), Matthews & Reid (2007) detected Hi emission co- 
incident in both position and velocity with RS Cnc, and 
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reported emission close to, but not unambiguously related 
to, EP Aqr. They also reported Hi emission from the cir- 
cumstellar shell of R Cas. 

Generally the H I emissions that we detect with the NRT 
are extended, indicating shell sizes ~ 1 pc. In some cases, e.g. 
EP Aqr (Paper I), they reveal complex spatial and dynamic 
structures. The 21 cm line thus probes regions much larger 
than those that can be studied with molecular lines, which 
are limited by photo-dissociation (r < 0.1 pc). Extended 
shells have also been observed in the far-infrared contin- 
uum emission (60 and 100 /j,m) by the Infrared Astronomical 
Satellite (IRAS, Young et al. 1993a). The spatial extension 
estimated from our Hi observations are comparable to, or 
possibly larger than, those obtained by Young et al. (1993a). 
Some of the infrared shells appeared detached from the cen- 
tral stars. The infrared source associated to Y CVn was 
reobserved by the Infrared Space Observatory (ISO), which 
clearly detected at 90 /im a detached dust shell surrounding 
the central star (Izumiura et al. 1996). These detached shells 
were interpreted as resulting from the interaction of the ex- 
pelled material with the Interstellar Medium (ISM, Young 
et al. 1993b) , or from a past event of mass loss much higher 
than today (Izumiura et al. 1996). However the far- infrared 
data lack crucial kinematic information. 

The H I line profiles obtained with the NRT have in gen- 
eral a gaussian shape, with a Full Width at Half-Maximum 
(FWHM) smaller than, or equal to, the width of the CO 
rotational lines, implying that the expansion velocity de- 
creases outward (Papers I & III). This slowing-down of the 
outflows indicates an interaction with local material, either 
of circumstellar or interstellar origin. These high-spectral- 
resolution H I profiles provided the first direct evidence of a 
slowing-down of circumstellar matter in the external shells 
of red giants. 

In the present paper, we revisit the case of Y CVn, 
because, as the ISO observations show images that are rather 
rounded, we expect that a spherical symmetry may apply 
and therefore that the interpretation will be easier. Also 
the source is relatively bright in H I, with a low interstellar 
contamination, so that in principle H I profiles of high quality 
are accessible. With these favorable conditions we expect 
that the H I profile modelling will bring useful constraints on 
the detached shell characteristics, and more generally on the 
detached shell phenomenon. We have acquired new Hi data 
on Y CVn in order to improve those already presented in 
Paper I. 



2 Y CVN 

2.1 The central star 

Y CVn is a J-type carbon star, i.e. a carbon star whose pho- 
tosphere is enriched in 13 C. The abundance ratio, 12 C/ 13 C 
~ 3.5 (Lambert et al. 1986), is close to the CNO-cycle equi- 
librium value (3-3.5, Lattanzio & Forestini 1999). The ef- 
fective temperature is ~ 2760 K (Bergeat et al. 2001), so 
that atomic hydrogen should dominate over molecular hy- 
drogen in the atmosphere and eventually above it (Glassgold 
& Huggins 1983). However, Lambert et al. (1986) detected 
a weak 1-0 S(0) H2 line at 2.22 /im, so that some hydrogen 
should still be in molecular form within the photosphere. 



The abundance analysis performed by Lambert et al. shows 
also that Y CVn may have a slightly sub-solar metallicity 
(— 0.2±0.2), with oxygen showing the most underabundance 
(-0.4). 

The exact stage of evolution of Y CVn is not clear. 
J-type stars do not show an enhancement in s-process ele- 
ments, suggesting that they have not gone through thermal 
pulses and that the surface carbon enrichment is not due to 
a third-dredge-up event. For instance, they might owe their 
carbon-rich composition to an He core flash when they were 
still on the Red Giant Branch (Dominy 1984). 

The distance derived from the Hipparcos parallax (4.59 
± 0.73 mas, Perryman et al. 1997) is ~ 218 pc. This places Y 
CVn at 207 pc from the Galactic Plane in a direction close to 
the Galactic North pole (b 11 = 72°). We adopt the Hipparcos 
distance, but note that it might be underestimated. Knapp 
et al. (2003) have re-reduced the Hipparcos data and give 
a revised parallax (3.68 ± 0.83 mas) that tends to place 
the source somewhat further from the Sun (272 pc). Also 
Bergeat et al. (2002) estimate the parallax at ~ 3.85 mas 
(or ~ 260 pc). 

The K magnitude is -0.74 (Guandalini et al. 2006) 
which translates to a bolometric magnitude of 1.96 (Le 
Bertre et al. 2001) or ~ 6200 L© at 218 pc. Such a luminos- 
ity places Y CVn clearly on the AGB, probably on the early 
AGB (E-AGB) rather than on the thermally pulsing AGB 
(TP- AGB). 

The proper motion observed by Hipparcos is -2.20, in 
Right Ascension (RA), and 13.05 mas yr~ , in Declination. 
When corrected for the solar motion towards apex, it trans- 
lates to 13.4 and 17.4 mas yr -1 respectively, and for a dis- 
tance of 218 pc to 14 and 18 kms -1 , with respect to the 
(solar) Local Standard of Rest (LSR). It corresponds to a 
motion, in the plane of the sky, towards North-East (posi- 
tion angle, PA = 38°). 

2.2 The circumstellar envelope 

Evidence for the presence of a circumstellar shell around 
Y CVn was found by Goebel et al. (1980) who detected the 
11.3-/im band ascribed to SiC dust grains. The 1-0 rotational 
line of CO was first detected by Knapp & Morris (1985). 

Knapp et al. (1998) have obtained a high-resolution CO 
spectrum of the wind from Y CVn. They determine a LSR 
radial velocity, Vi sr = 21.1 kms -1 , an expansion velocity, 
V CX p= 7.8 kms -1 , and a mass loss rate, M = 1.1 10~ 7 
M Q yr -1 (at 218 pc). Schoier et al. (2002) obtained sim- 
ilar estimates with V eX p = 8.5 kms -1 and M = 1.5 10 -7 
M0yr _1 , as well as Teyssier et al. (2006) with V cxp = 6.5 
kms -1 and M = 1.0 10~ 7 M Q yr -1 . The CO emission is ex- 
tended with diameters, (j> ~ 13" in the 1-0 transition, and 
~ 9" in 2-1 (Neri et al. 1998). A faint asymmetry is sus- 
pected in the 2-1 data. 

An intense emission in the 13 CO 1-0 line, detected by 
Jura et al. (1988), confirms the large abundance of 13 C in 
the Y CVn outflow. Schoier & Olofsson (2000) evaluate the 
12 CO/ 13 CO abundance ratio to ~ 2.5, a ratio larger than the 
equilibrium value, indicating a possible enrichment in 13 C 
of the CO molecule within the outflow. The mass loss rates 
estimated from 1 CO rotational line data may therefore be 
underestimated by a factor ~ 1.4. We note also that the 
underabundance in oxygen, as suggested by Lambert et al. 
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Figure 1. Map of the 21 cm Hi emission from the Y CVn circumstcllar envelope observed with the NRT. The steps are 4' in RA 
(1 beam) and 11' in Declination (1/2 beam). North is up and East to the left. 



(1986), would directly translate into an underabundance in 
CO. 

An extended emission at 60 and 100 /im was discovered 
by IRAS (Young et al. 1993a). This emission could be mod- 
elled with a resolved isothermal shell, at ~ 30-40 K, of inner 
radius, ri„~ 2.8', and outer radius, r out ~ 5.5'. An image at 
90 (im with a better spatial resolution (~ 40") and a bet- 
ter sensitivity was obtained by ISO (Izumiura et al. 1996). 
It clearly shows a central infrared source surrounded by a 
quasi-circular detached shell, with r 4n ~ 2.8' and r ou t~ 5.1'. 
This image of a detached shell supports the modelling of 
the IRAS data by Young et al. (1993a). In addition, in the 
ISO image at 90 /im, the detached shell seems displaced 
to the West by about 0.5-1' with respect to the central in- 
frared source. Nevertheless, the circular morphology of the 
detached shell and the small size of the offset, as compared to 
the source diameter, leave no doubt on the association with 
Y CVn. Adopting a dust-to-gas mass ratio of 4.5 10~ 3 , the 
mass in the detached shell can be estimated at 0.06 M (at 
218 pc). 



This detached shell, observed in the dust continuum 
emission, is interpreted by Young et al. (1993b) as a prod- 
uct of the slowing-down of an expanding circumstellar shell 
by the surrounding ISM. On the other hand, for Izumiura 
ct al. (1996), it is the result of an intense episode of mass 
loss, about two orders of magnitude larger than presently 
(as probed by the CO emission, which is limited to the 
central ~ 10", Neri et al. 1998), with an outflow velocity 
~ 15kms _1 , that lasted about 2 10 4 years and stopped 1.4 
10 4 years ago. 

In Paper I, we reported the detection of Y CVn in the 
Hi line at 21 cm. The emission is extended compared to the 
NRT beam size in Right Ascension (4'). The line profile is 
composite with a spectrally narrow component (Comp. 1, 
FWHM ~ 2.9 kms -1 ) superposed on a weak broad com- 
ponent (Comp. 2, FWHM ~ 14 kms -1 ). The 2 compo- 
nents are centered at Vi sr ~ 20-21 kms -1 , in good agreement 
with the CO lines. Comp. 1 was found to be spatially ex- 
tended (~ 12') and to be more intense West than East. The 
H I source traced by Comp. 1 was associated to the detached 
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shell observed by ISO, whereas the other one was associated 
with the Y CVn outflow already observed in CO molecular 
lines. 

The spatially extended emission associated to Comp. 1 
and its quasi-gaussian line-profile of width ~ 3 kms -1 were 
interpreted as an indication that the outflow velocity should 
decrease outward, a result which agrees with the Young et al. 
(1993b) interpretation of a slowing-down by the surrounding 
ISM. However it does not exclude a brief episode of a large 
mass loss rate in the past (see Sect. [6]). 



3 OBSERVATIONS 

The NRT has a rectangular aperture of effective dimensions 
160 mx 30 m. At 21 cm, its beam has a FWHM of 4' in Right 
Ascension (RA) and 22' in Declination. The data on Y CVn 
have been acquired in the position-switch mode with off- 
positions taken at ±4', ±8', ±12', ±16', and ±32' from the 
central position in the East- West direction. The central posi- 
tions were selected on the position of the star (adopted from 
the SIMBAD database: 12 h 45 m 07.8 s ±45°26'25", 2000.0), 
and at ±11' and ±22' from it, in the North-South direction. 
The spectral resolution corresponds to 0.32 kms -1 . For more 
details we refer to Paper I. 

In order to obtain the best Hi line profiles needed for 
our modelling, we have re-processed our previous data with 
an improved procedure (better weighting of the individual 
observations) and acquired new observations, with a band- 
width twice as large, allowing a better determination of the 
baselines. In total, 107 hours of data have been acquired be- 
tween Sept. 2002 and Jan. 2007, and are used in the present 
work. 

The new map, that results from the merging of old and 
new data (Fig. [TJ, has a much better quality than the one 
presented in Paper I. The signal-to-noise ratio is on average 
twice better than in our 2004 map (Paper I). 

The H I line at the stellar position (full line in Fig. [2} 
is clearly dominated by a narrow quasi-gaussian component 
(Comp. 1) of FWHM ~ 3.1 kms -1 , and intensity ~ 360 mjy, 
superposed on a weak rectangular pedestal of width ~ 16 
kms -1 , and intensity ~ 10 mjy. In Table Q] we give the 
results of the line fitting by the sum of a gaussian and of 
a rectangle; as this fitting closely matches the data, it is 
not represented in the figures. The dashed lines result from 
the modelling of the whole map that will be discussed in 
Sect. [5] The radial velocities of the two components are 
nearly the same and in agreement, within 0.5 kms -1 , with 
that of Y CVn as determined from CO data (Knapp et al. 
1998). The improved quality of our Hi spectrum allows to 
see that Comp. 2 might have a rectangular profile rather 
than a gaussian one (cf. Paper I); nevertheless, for continu- 
ity, we keep the same notation, Comp. 1 & 2, as in Paper 
I. 

East and West of the central position (Fig. [3] upper 
panel) the Hi line profile is essentially gaussian-like, i.e. with 
the present signal-to-noise ratio we do not detect a pedestal. 
We confirm the asymmetry between East and West, which is 
consistent with that seen on the ISO image at 90 /im. From 
the relative fluxes of Comp. 1 between the central position 
and the two offset positions, we estimate the displacement 
at ~ 1' West. As noted in Sect. 12. ll this displacement is small 
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Figure 2. Hi spectrum obtained with the NRT on the stellar 
position (thin line) and modelled spectrum (Sect. [5] dashed line). 
A zoom on the lower part of the spectrum is displayed in the 
lower panel in order to show better the pedestal (Sect. [3). The 
vertical dotted line indicates the radial velocity adopted for the 
model (Visr= 20.6 kms -1 ). 



Table 1 . Fits to the two H I components of Y CVn (not repre- 
sented in the figures). Comp. 1 is fitted with a gaussian profile, 
and Comp. 2 with a rectangular profile. In the case of Comp. 1 we 
quote the FWHM and for Comp 2 the full width. For comparison 
the CO (3-2) profile obtained by Knapp et al. (1998) gives Vi sr = 
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compared to the source size. In addition, we note that the 
radial velocity of Comp. 1 is consistent with that of Y CVn 
further supporting their association. We do not detect the 
source at ±8' from the stellar position (Fig. [3] lower panel). 
The diameter of the Hi source can then be estimated at 
8' ±4'. 

The source is clearly detected at ±11' in the North- 
South direction (Fig. [4] upper panel) with more than half the 
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Radial velocity (km/s) 

Figure 3. Upper panel: average Hi spectrum obtained with the 
NRT on the two off-positions at ±4' in RA (thin line) and mod- 
elled spectrum (Sect. [5] dashed line). Lower panel: same, with the 
two off-positions at ±8'. 
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Figure 4. Upper panel: average H I spectrum obtained with the 
NRT on the two positions at ±11' in Declination (thin line) and 
modelled spectrum (Sect. [5] dashed line). Lower panel: same, with 
the four off-positions at ±11' in Declination and ±4' in RA. 



flux on Y CVn. Therefore it is also extended in Declination, 
but our spatial resolution is not sufficient to characterize the 
size in this direction. In the map, we note that the source is 
slightly more intense North than South, by about 10 % at 
±11', which is an indication of a slight (relative to the beam) 
displacement North that we estimate at l'±0.5'. Finally, the 
source is also detected at the four positions ±4' in RA and 
±11' in Declination (Fig. 3] lower panel). 

By integrating the line profiles over the full map, the 
total Hi flux associated to Y CVn can be estimated at 3.38 
Jykms" 1 . At 218 pc, this translates to 3.8 ± 0.2 10" 2 M Q in 
atomic hydrogen. Most of it is associated to Comp. 1. We 
estimate the flux corresponding to Comp. 2, from the width 
and height of the pedestal detected on the Y CVn position, 
at 0.16 Jykms" 1 , or 1.8 10" 3 M . 



4 INTERPRETATION 

Although we are limited by the spatial resolution of the NRT 
(4' x 22'), by combining our high-spectral resolution H I data 
with the ISO image obtained in the dust continuum by Izu- 
miura et al. (1996), we can obtain a description of the Y 
CVn circumstellar environment that takes into account the 
kinematics of the gas. 

The source traced by Comp. 2 can be associated with 



the central infrared source. Its rectangular profile is typi- 
cal of an optically thin H I emission of a uniformly expand- 
ing shell (Paper I). The expansion velocity (~ 8 kms -1 ) is 
consistent with the outflow velocity obtained from CO ro- 
tational lines. As Comp. 2 is also centered on the radial ve- 
locity of Y CVn, we can infer that the corresponding source 
has a nearly spherical geometry. 

The source traced by Comp. 1 is spatially resolved by 
the NRT and has a diameter, (f> ~ 8' ± 4'. Its hydrogen mass 
is consistent with the total mass of the detached shell esti- 
mated by Izumiura et al. (1996). Therefore it is reasonable 
to associate it to the infrared detached shell. The spectral 
width, ~ 3 kms -1 , shows that the material is moving out- 
wards at a reduced velocity as compared to the source traced 
by Comp. 2. This is in agreement with the analysis of Young 
et al. (1993b). 

We now pursue this analysis and consider the interac- 
tion of the Y CVn outflow with external matter. The in- 
frared image and our H I map show an offset of 1' West with 
respect to the central source, but this offset is small com- 
pared to the diameter of the detached shell, and in the ISO 
image the detached shell has a quasi-circular appearance. 
Also the central velocity in Hi (20.6 kms -1 ) is close to that 
in CO (21.1 km s -1 ), so that no large deviation along the line 
of sight can be detected. In the following we assume spheri- 
cal symmetry, which in particular implies that the external 
medium is at rest with respect to Y CVn. 
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Figure 5. Schematic view of the Y CVn detached shell. The 
termination shock is located in n, the contact discontinuity in 
Tj, and the bow shock in r<z. CS stands for circumstellar material 
and EX, for external material (see Sect. [I). 



The interaction of a spherical outflow with external 
matter leads to the formation of a region of compressed ma- 
terial within two spherical boundaries (Fig. [5] Lamers & 
Cassinelli 1999). The internal boundary, n, defines the sur- 
face at which the supersonic outflow is abruptly slowed-down 
by compressed circumstellar material (termination shock). 
The external boundary, r 2 , defines the surface at which the 
external medium is compressed by the expanding shell (bow 
shock). Within these two limits we find compressed materi- 
als of stellar (CS) and external (EX) origins that are sepa- 
rated by a contact surface defined by rf. In this description 
no material is allowed to flow over the contact discontinuity 
in rf. Finally, inside r\, the stellar outflow is in free expan- 
sion. 

We interpret the detached shell observed by ISO and 
the Hi source traced by Comp. 1 as this region of com- 
pressed material delineated by r\ and r 2 , and identify n 
with Tin and r 2 with r ott t. I.e. we assume that the dust and 
the gas have a similar spatial distribution; this assumption 
is further discussed in Sect. [6] We also adopt the best esti- 
mates from Izumiura et al. (1996), so that for a distance of 
218 pc, ri = 0.18 pc, and r 2 = 0.32 pc. 

The width of the H I emission allows to estimate an up- 
per limit to the average Hi temperature within the com- 
pressed shell. For a maxwellian distribution of hydrogen 
atoms at temperature T, the Hi emission projected on a 
line of sight has a gaussian profile of width (Gerard 1990): 



FWHM(km s" 1 ) = 0.214(T(iT)) 1/2 



(1) 



Therefore, the temperature of the gas within the detached 
shell is, on average, at most 210 K. Conversely, the bulk of 
the material in the detached shell is moving outwards at a 
maximum velocity ~ 1/2 FWHM, i.e. ~ 2 kms -1 . 

The source associated to Comp. 2 can be identified with 
the freely flowing stellar wind that fills the region inside the 
termination shock (r < r%). The crossing time is about 22 
10 3 years, which translates to an average mass loss rate in 
atomic hydrogen of about 0.8 10~ 7 Mq yr _1 . 

The atomic hydrogen mass of the detached shell, Mjjt 
= 3.6 10~ 2 M S) 

can be separated in two components : a cir- 
cumstellar part, Mdt,cs, located between n and rf, and an 
external one, Mdt,ex, between rf and r 2 (see Fig. [5}. To es- 



timate the latter we assume that it is equal to the hydrogen 
ISM mass enclosed in a sphere of radius r%. Following Young 
et al. (1993b) we adopt an ISM density of 0.25 H cm -3 from 
a galactic scale height of 100 pc, and an ISM density of 2 H 
cm -3 in the galactic plane. In these conditions, Mdt,ex ~ 
0.9 10 -3 Mq, and therefore M_dt,cs can be estimated at ~ 
3.5 10 -2 Mq in atomic hydrogen. The mass of the detached 
shell appears clearly dominated by circumstellar material. 
We note that Y CVn is relatively close to the Sun, in a di- 
rection for which the interstellar medium might be deficient 
(Sfeir et al. 1999), so that we might have been overestimat- 
ing the ISM density. 

Finally, if we assume that the mass loss rate has been 
constant with time, and adopt the total mass derived from 
Comp. 2, Mdt,cs, we obtain a characteristic time for the 
formation of the detached shell of tos ~ 4.5 10 years, ~ 
20 times larger than the crossing time of the freely flowing 
wind region. A characteristic velocity for the expansion of 
the detached shell is then: ri/tos ~ 0.4 kms -1 . This is 
probably an overestimate of the present expansion rate, as 
the expansion rate should decrease with time (Young et al. 
1993b). 

In conclusion the detached shell appears as a region 
of the Y CVn environment where circumstellar matter has 
been slowed down and is stored for a few 10 5 years before 
being injected in the ISM. 



5 MODELLING 

For the freely expanding wind zone (r < ri), we consider an 
outflow with a constant velocity (uq = 8 kms -1 ) and a vo- 
lumic density varying in 1/r 2 , corresponding to a constant 
hydrogen mass loss rate. The temperature of the gas is as- 
sumed to be low enough that its effect on the H I profile can 
be neglected. Indeed, at distances larger than 10 16 cm the 
gas temperature should be lower than 100 K (e.g. Kahane 

6 Jura 1994) and the thermal broadening should stay small 
compared to vo. 

In order to describe the detached shell around Y CVn 
(ri < r < r 2 ) in simple terms, we make the hypotheses of 
stationarity and spherical symmetry; all quantities depend 
only on r, the distance to the central star, and in particular 
the limits n, r/ and r 2 are treated as independent of time. 
The mass flow is constant, and the same rate apply in the 
detached shell from ri to r / as in the freely expanding wind 
zone (r < ri): 



M = 47rr pv 



(2) 



where p and v are the density and velocity in r. For M , we 
adopt the value corresponding to that determined in Sect.|4j 
Mm = 0.8 10 -7 Msyr 1 . 

We assume that the gas behaves as an ideal gas of mean 
molecular weight, /j,. The pressure (p), density and temper- 
ature (T) are then related by the equation of state: 

kT 2 ,„\ 

P = P = pc (3) 

pm H 

where c is the isothermal sound velocity, k the Boltzmann 
constant, and tuh the mass of the hydrogen atom. For a 
neutral atomic gas with 10 per cent 4 He, and 90 per cent H, 
H = 1.3. 



The detached shell around Y CVn 7 



The equation of motion for an ideal gas in spherical 
geometry is given by : 

dv 1 dp Ik dpT . 
dr p dr p prriH dr 

(i.e. we neglect the stellar gravity). Using Equ. [2] which is 
valid from ri to r/, we re- write Q as: 



dr v r 



k dT 



(5) 

firriH dr 

The properties of the gas in n are given by the jump 
conditions (Dyson & Williams 1997): 



povo = pivi 
vi Mi + 

('() 



(6) 
(7) 



where po, vq and Mo are the upstream density, velocity and 
Mach number. Equ. 6 is a general expression while Equ. 7 
only applies to a mono-atomic adiabatic gas. From the pre- 
vious discussion we adopt vo ~ 8 kms -1 and deduce po 
from the mass loss rate. We estimate Mo ~ 20 assuming an 
upstream temperature, To — 20 K. The downstream tem- 
perature within ri, Ti, is given by: 



3fim H 2 



isook 



(8) 



As the average temperature within the detached shell is at 
most 210 K (previous section), the gas must cool down to 
a temperature lower than 210 K within the time lapse of 
the detached shell formation (~ 4 10 s years). This cooling 
should occur through atomic lines and/or dust emission, and 
probably depends on metallicity. (Some heating of the gas by 
dust may also be expected close to ri, see Sect. [6]) We do not 
estimate the cooling rate, but rather adopt a temperature 
profile through the detached shell that we will constrain from 
our Hi observations. In practice we adopt a law of the form: 



T r 
log — =a log — 

with Ti = 1800 K. Equ. ([SJ can then be re-written as: 



dv 
dr 



(v 



c 

V 



(9) 



(10) 



The equation of motion (|10p is solved layer by layer, 
using the classical Runge-Kutta method, outwards from n, 
under the initial conditions of density, temperature and ve- 
locity, pi, Ti and Vi. The limit rj is set such that the hy- 
drogen mass between n and r/ is equal to Mdt,cs- 

From rf to r2, we cannot use Equ. [5] As the amount 
of interstellar matter between r/ and r2 is small compared 
to the mass in the detached shell, we simply assume a 1/r 2 
dependence for the density in this region, with the condition 
that its hydrogen mass is equal to Mdt,ex, i.e. 0.9 10~ 3 
Mq. The temperature is assumed constant from if to r2 
(although strictly speaking it is expected to increase from 
if to T2)- The velocity is then obtained from Equ. [4] 

The velocity, hydrogen-density and temperature profiles 
are then used in an Hi emission model that we have al- 
ready developed (Papers I and II). For each point within the 
detached shell, the thermal broadening is derived through 
Equ. [1] In the model the Hi emission from an envelope of 
matter flowing radially from the central star is convolved 




2 4 
radius (arcmin) 



radius (arcmin) 

Figure 6. Upper panel: adopted temperature profile in the de- 
tached shell for the model described in Sect. [5] Lower panel: ve- 
locity profile. The dashed line represents the isothermal sound 
velocity (c). 



with the response given by a telescope of rectangular aper- 
ture with effective dimensions 160 mx 30 m. We thus calcu- 
late the H I emission as it should be observed by the NRT at 
the central position and at different positions offset in RA 
and Declination (dashed lines in Fig. [2] [3] and |4)| . It should 
be noted that in the present modelling the temperature pro- 
file plays a double role, first in the kinematics through the 
isothermal sound velocity (Equ. [3j , second in the H I profile 
through thermal Doppler broadening. 

Our modelling depends only on the value selected for 
the constant a in Equ. [9] For the comparison with the obser- 
vations (Fig. [2] and following) we adopted a stellar velocity, 
Vi sr = 20.6 kms -1 . We tested several values and obtained a 
satisfactory fit to the observed Hi line profiles for a = —6.0. 
The resulting temperature profile in the detached shell is 
plotted in Fig. [6] (upper panel); it reaches a minimum in r/ 
at ~ 165 K. The corresponding velocity profile is shown in 
the same figure (lower panel) . Within the detached shell, the 
flow remains sub-sonic; the velocity profile that we obtain 
is similar to the "breeze" solution of the classical theory of 
stellar winds (e.g. Lamers & Cassinelli 1999). The density 
profile and the flux of matter are shown in Fig. [7] Both are 
clearly peaked at 77 = 0.27 pc (4.22'). 

The numerical values used for the parameters of the 
model are summarized in Table [2] 

The Hi emission model developed in Papers I and II 
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radius (arcmin) 




radius (arcmin) 

Figure 7. Upper panel: atomic hydrogen density profile for the 
model described in Sect. [5] Lower panel: atomic hydrogen mass- 
flow profile. The vertical dotted lines mark the radii, r\, rj and 
V2, of the model. 



assumes that the brightness temperature is proportional to 
the atomic hydrogen column density (i.e. T ^ 10 K) and 
that the emission remains optically thin (r <C 1). The first 
hypothesis is clearly verified (see Fig. H3 upper panel). The 
optical depth of the Hi line is given by: 



3c 



light 



1 „ »r hv 1 

-^AwNb — ^— 



32tt i/ 2 'J:TAv 
where cu g ht is the velocity of light, Am the spontaneous 
emission coefficient (2.8688 10 -15 s -1 ), N H the atomic hy- 
drogen column density, h the Planck constant and Av the 
line width. Expressing the line width in kms -1 (AV) and 
the column density in cm -2 : 



r = 5.50 10" 



TAV (12) 

From Fig. [8] the column density is maximum for r = 3.72' 
(0.24 pc). Since T > 160 K and AV ~ 3 kms -1 , the optical 
depth is always smaller than 0.03 across the detached shell. 



6 DISCUSSION 

Our modelling, that assumes stationarity and sphericity, is 
simplified. However, it allows to grasp the physical condi- 
tions within the Y CVn detached shell and already pro- 
vides a fair adjustement to the high-resolution Hi spectra 




radius (arcmin) 

Figure 8. Atomic hydrogen column density in the detached shell 
model. 



Table 2. Model parameters (d = 218 pc). 



M (in hydrogen) 


0.78 10" 7 Meyr 1 




1.3 


ti 


22 500 years 


tus 


427 500 years 


ri 


0.18 pc (2.84') 


T f 


0.27 pc (4.22') 




0.32 pc (5.10') 


T (= T"), T+ 


20 K, 1804 K 


T/ (= T a ) 


167 K 


vo(s v 1 ), v+ 


7.8 kms" 1 , 1.96 kms 


Vf 


0.065 kms" 1 


V2 


0.9 kms" 1 


n", n+ 


1.0 Hem- 3 , 3.9 Hem 


n f , n+ 


53. Hem" 3 , 0.7 Hem 


ri2 


0.5 Hem" 3 


M r < ri (in hydrogen) 


1.75 10" 3 M 


M DT,CS ( in hydrogen) 


3.32 10- 2 M 


Mot, ex (in hydrogen) 


0.87 10" 3 M 



obtained at different positions with respect to the central 
star. The fits to some line profiles appear too broad, in par- 
ticular on the central position (Fig. [3]) for which we have 
the highest signal to noise ratio. However, the total Hi flux 
in the map is, by construction, well reproduced. It could be 
possible to improve the fits by decreasing the mass loss rate 
in the model, but this would be at the expense of the agree- 
ment with the mass loss rate determined from CO rotational 
lines, and with the total H I mass. 

In Fig. [2] we note a shift in velocity of ~ 1 kms -1 be- 
tween the model pedestal and the observed one around 29 
kms - . It comes from the choice in the model of a central 
velocity at 20.6 kms -1 , whereas the stellar velocity is more 
likely at 21.1 kms -1 , as can be estimated from the CO radial 
velocity (e.g. Knapp et al. 1998). As Comp. 1 is well fitted 
in velocity, it indicates that the detached shell is on average 
moving at a velocity of ~ 0.5 kms -1 with respect to the 
central star. This could be an effect of a systematic motion 
of the star with respect to surrounding matter that would 
exert an asymmetrical pressure on the stellar outflow. This 
effect, observed along the line of sight, could be related to 
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the East- West asymmetry observed in our H I map (Fig. [TJ 
and in the ISO image at 90 (im (Izumiura et al. 1996). We 
note that the Y CVn proper motion towards East (PA = 
38°, Sect. I2.1|l is consistent with this East- West asymme- 
try. On the other hand we have presently an indication for 
a North-South asymmetry that should be investigated with 
better spatial resolution. 

In order to test the assumption that the detached shell 
of Y CVn is displaced by 1' West, we have performed a 
new modelling of the Hi profiles assuming that the source 
is 1' West of the centre of our observation grid, all other 
model parameters (Table [2} being kept identical. The re- 
sults are shown in Fig. [§] The intensities 4' East and West 
of Y CVn are well reproduced and the agreement on the cen- 
tral position is even improved. This is an indication that the 
displacement is real. In that case the hypothesis of spheric- 
ity should be re-examined and a non-spherical model (that 
might also account for the 0.5 kms -1 velocity shift discussed 
above) should be considered. 

A surprising feature of our model is the relatively low 
density between rj and r-z that is not confirmed by the ISO 
image. As we do not detect any H I emission at 8' East and 
West of Y CVn (Fig. [3] lower panel), we suspect that the 
outer radius of the H I detached shell has been overestimated 
by our hypothesis that hydrogen and dust have the same 
spatial distribution, i.e. that ri — r out — 5.1'. The displace- 
ment of the infrared source West of Y CVn may have led 
also to an overestimate of r ou t by Izumiura et al. (1996). 

In our model we adopted an empirical temperature pro- 
file. Ideally it would have been better to determine the tem- 
perature from a cooling function. This would have substan- 
tially increased the complexity of our modelling because the 
temperature enters into the equation of motion (Equ. [5}. 
Also the cooling function is expected to depend on metal- 
licity which is not well constrained for Y CVn. In addition 
the mean free path of dust particles through the gas should 
be ~ 10 9 cm, so that some heating of the gas due to friction 
should occur, especially close to n where the relative veloc- 
ity is ~ 6 kms -1 . The grain photoelectric effect may also 
contribute to heating. The adoption of an empirical profile 
for the temperature includes all these effects. We note that 
a shallower profile would have led to a too broad Hi emis- 
sion profile, and that, with a steeper one, the model would 
have failed to reproduce a correct emission on the East /West 
position. 

The hydrogen mass loss rate that is deduced from the 
pedestal and from the inner radius of the detached shell 
is ~ 0.8 10 -7 M0yr -1 . It translates to a total mass loss 
rate of ~ 1.0 10 -7 Mgyr -1 , a value compatible with those 
derived from CO line fitting. However, as noted in Sect. 12.21 
the latters are probably underestimated by a factor 1.4 or 
more. Considering the uncertainty on the CO abundance, a 
mass loss rate derived from Hi should in principle be more 
reliable. However we note that our estimate corresponds to 
an average over 22 10 3 years or even longer. It does not 
exclude possible variations on smaller timescales. On the 
other hand the CO estimates correspond to the last period 
of ~ 10 3 years. From a detailed CO modelling Dinh-V-Trung 
& Nguyen-Q-Rieu (2000) concluded that the Y CVn mass 
loss rate must have increased by a factor 2 or more over the 
last 1600 years. Y CVn may thus be presently in an episode 



of high mass loss as compared to the average value that we 
obtained in Hi . 

In our model of the detached shell we assume that the 
mass loss rate has been approximately constant for about 
4-5 10 years. This may seem too large a duration in view 
of the expected time between two thermal pulses (~ 10 5 
years or less for carbon stars, Straniero et al. 1997) . However 
we noted that Y CVn is a carbon star of the rare J-type 
and that it may not necessarily be on the TP-AGB. On 
the other hand the properties that we implicitly assume for 
the Y CVn outflow fit rather well with the B-type models 
calculated by Winters et al. (2000) for carbon-rich late-type 
giants with stellar temperature ~ 3000 K and low variability. 
These models develop stationary outflow structures with no 
strong temporal variations. 

Our H I data modelling confirms that the detached shell 
of Y CVn is a dense (~ 50 cm -3 ), relatively mild and quiet 
(~ 200 K, ~ 1-2 kms -1 ) region where circumstellar matter 
stays for at least 4 10 years before being mixed with the 
ISM. 

Recently, using far-infrared data from the Spitzer Space 
Telescope, Ueta et al. (2006) have detected a bow shock neb- 
ula with a cometary shape around the AGB star R Hya. 
Wareing et al. (2006) have modelled this nebula in terms of 
an interaction of the stellar wind with the surrounding ISM 
through which R Hya is moving. They estimate the temper- 
ature behind the shock at ~ 35 000 K, which corresponds to 
a star-ISM relative velocity of ~ 35 kms -1 . At such tem- 
perature, hydrogen is expected to be entirely ionized. 

For Y CVn, although the transverse motions and the 
radial velocity indicate a motion relative to the ISM of com- 
parable magnitude (~ 31 kms -1 ), we find only a weak evi- 
dence for a distortion that could be ascribed to this motion. 
We cannot exclude the presence of ionized hydrogen between 
the contact discontinuity, r/, and the bow shock, ri. How- 
ever, we note that the mass in atomic hydrogen that we 
estimate from our Hi data, 0.04 Mg, agrees with the total 
mass estimated from the dust continuum by Izumiura et al. 
(1996), ~ 0.06 M . Therefore the detached shell around Y 
CVn seems to be dominated by compressed neutral matter. 
Our modelling assumes spherical symmetry, i.e. that the sur- 
rounding material is at rest with respect to the central star. 
This suggests that the external matter, within which the de- 
tached shell is expanding, is moving together with Y CVn. 
This material might therefore not be genuine ISM, but pos- 
sibly material remaining from an older episode of mass loss. 
We recall that the stellar evolution models (e.g. Renzini & 
Fusi Pecci 1988) invoke an important mass loss, perhaps ~ 
0.2 Mq, on the Red Giant Branch (RGB), that has not been 
detected up to now. In such a case our way of estimating the 
mass in the swept-up region (from r/ to r-i) % Mdt,ex ~ 0.9 
10 -3 Mq in atomic hydrogen (Sect. [3]), which relies on the 
assumption that external matter is interstellar, should be 
revised. For instance an increase of Mdt,ex by a factor 10 
(~ 0.9 10 -2 Mq) would imply a reduction of Mdt,cs to ~ 
2.5 10 -2 Mq and a reduction of tns to ~ 3.3 10 5 years. This 
would somewhat change the structure of the detached shell 
in the model, but not fundamentally. 

The phenomenon observed in Y CVn seems different 
from that leading to the formation of the detached CO shells 
observed around some carbon stars (Schoier et al. 2005). 
These molecular shells have a typical radius ^ 0.1 pc and 
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Figure 9. Effect of a shift 1' West of the Hi source with respect to the centre of the observation grid. Left panel: Hi spectrum obtained 
4' East; centre panel: Hi spectrum obtained on the centre; right panel: Hi spectrum obtained 4' West. The modelled spectra (dashed 
lines) are discussed in Sect. \E\ 



are expanding at a velocity ^ 12 kms" 1 (e.g. Olofsson et 
al. 2000). They are smaller than the Y CVn detached shell 
(also they look thinner with a CO diameter-to-thickness ra- 
tio ~ 5-10) and do not seem to be slowed down efficiently by 
an external medium, although there are signs of interaction 
with a surrounding medium. The sources with CO detached 
shells show complex CO line profiles with a double-peak (e.g. 
Olofsson et al. 1990) which is not seen in any of the Y CVn 
CO spectra (Knapp et al. 1998; Teyssier et al. 2006). Olof- 
sson et al. (1990) have suggested that these CO detached 
shells are produced by a mass loss eruption, possibly initi- 
ated by an He-shell flash. This mecanism is unlikely to play 
a role for Y CVn, because the central star is probably still 
on the E-AGB. 

The dynamical evolution of the interaction between 
stellar winds and the ISM has been explored by Villaver et 
al. (2002) . They obtain transient shells (~ 20 000 years) as- 
sociated with the wind variations induced by thermal pulses. 
Longer-lived shells could occur via interaction between suc- 
cessive events of enhanced mass loss, or via the continu- 
ous accumulation of ejected material in the interaction re- 
gion with the ISM. Mattson et al. (2007) have modelled the 
formation of detached shells around TP-AGB stars. They 
need to combine an eruption of mass loss to the interaction 
with external matter in order to explain the properties of 
CO detached shells. They obtain thin detached shells with 
small relative thickness (AR/R ^ 0.01). From their mod- 
elling, the cases for which the observed relative thickness is 
large (AR/R ^ 0.1) may have another origin. Therefore the 
mecanism for the formation of the Y CVn detached shell 
might be different from that responsible for the molecular 
detached shells around carbon stars. Given our low spatial 
resolution we cannot exclude short time (~ 100 years) fluc- 
tuations of the Y CVn mass loss rate such as those found, 
for instance, for the carbon star IRC +10216 (Mauron & 
Huggins 2000). 

On the other hand, if the scenario that we propose for 
Y CVn is correct, we expect to find similar detached shells 
around other AGB stars, and in particular around mass- 
losing oxygen-rich ones, because there is no reason for it to 
be specific to carbon stars. From our survey of atomic hy- 
drogen around evolved stars (Paper III) we find H I emission 
from oxygen-rich as well as from carbon-rich stars. Also the 
H I line profiles that are observed indicate a slowing-down of 
circumstellar matter at large distance from the central stars. 



7 PROSPECTS 

The spatial resolution of our Hi data is limited to ~ 4' in 
RA. Therefore we have used the ISO image at 90 /xm to 
better constrain the geometry of the circumstellar shell. On 
the other hand the high spectral resolution of our H I data 
has allowed to study directly the kinematics and the physical 
conditions within the Y CVn detached shell. Also the large 
collecting area of the NRT has been useful to reach a high 
sensitivity on the low surface brightness Hi emission of Y 
CVn. 

However our interpretation remains limited by the low 
spatial resolution in Hi . For instance, we have used the in- 
frared image as a guide and identified n (Hi) with Ti n (dust) 
and T2 with x ou t. It implicitly assumes that the dust and the 
gas spatially coincide. Also the infrared emission depends 
on the dust temperature which should vary with distance to 
the central star. An imaging in H I would allow to constrain 
directly the gas density within the shell as a function of the 
distance to the central star. Indeed in our model we have 
shown that the optical depth is always smaller than 1 and 
that the temperature in the detached shell is larger than 
10 K, so that the Hi emission at each position should be 
directly proportional to the column density. Therefore in- 
terferometric data in Hi with a high sensitivity and a high 
spectral resolution would be very valuable to probe the gas 
physical conditions and kinematics within the detached shell 
of Y CVn. 

Such high angular resolution imaging in H I should allow 
to determine clearly n and to compare it to r; n . It should 
also allow to determine the offset of the H I detached shell 
with respect to Y CVn. The spectrally broad component (2) 
is expected to be confined to the interior of the gas detached 
shell, which is expected to be seen only in the spectral band- 
width corresponding to the narrow component (1). Finally, 
we predict a concentration of hydrogen at r ~ 0.24 pc (3.7', 
see Fig. |8j). 

Other atomic lines, for instance C I at 492 GHz, might 
also be valuable tracers of the gas within the detached shell. 
Cl was not detected by Knapp et al. (2000); however, they 
observed on the stellar position with a beam of only 15" and 
a throw of 60" (< ri). In the optical range, resonant lines 
of Na or K could be useful. Mauron & Guilain (1995) did 
not detect any Nal/Kl emission, but similarly they looked 
only close to the central star (5"). 

The conditions for the temperature and the density that 
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we find in our model might be favorable to the formation of 
molecular species through a non-equilibrium chemistry as it 
has been suggested for some regions of the diffuse ISM (Fal- 
garonc ct al. 2005) . In fact as the physical conditions within 
the detached shell and the timescale can be characterized 
through modelling, Y CVn could thus be a good target for 
studying this type of ISM chemistry. 



comments that helped us to improve this paper. We thank 
Jean Borsenberger for developing a new and efficient pro- 
cedure for processing NRT data, and Jerome Pety for his 
advices in using the GILD AS environment. This research 
has made use of the SIMBAD database, operated at CDS, 
Strasbourg, France and of the NASA's Astrophysics Data 
System. 



8 CONCLUSIONS 

We have presented high spectral resolution H I data obtained 
on Y CVn with the NRT. The emission is spatially resolved 
with a diameter ~ 8' ± 4'. The spectrum obtained on the 
stellar position reveals a rectangular pedestal centered at 
+21.1 kms -1 (Comp. 2) that traces an 8 kms -1 outflow of 
~ 1.0 10~ 7 M Q yr _1 . This outflow corresponds fairly well to 
the wind detected in the CO rotational lines. 

The spectrally narrow component (Comp. 1) is centered 
at +20.6 kms -1 and traces an expanding gas shell that we 
associate with the detached shell imaged by ISO in dust 
continuum emission at 90 /im (Izumiura et al. 1996). The 
bulk gas temperature within this shell is ~ 100-200 K and 
the expansion velocity ~ 1-2 kms -1 . We have developed a 
simplified model in which the detached shell is produced by 
a slowing-down of the Y CVn outflow by external matter. In 
this model the mass loss is taken constant for ~ 5 10 5 years 
with the same rate as presently, i.e. 1.0 10~ 7 Mq yr _1 . The 
sharp decrease of the velocity by a factor ~ 4 is interpreted 
as due to a shock at the inner boundary (termination shock). 
The Y CVn detached shell thus appears to act as a lock 
chamber where circumstellar matter is stored for a few 10 s 
years before being injected in the ISM. 

We note that in our Hi model of the detached shell, 
we have adopted the same dimensions as those of the dust 
detached shell. These dimensions might not apply exactly 
and a high angular resolution mapping at 21 cm is needed. 
It would also provide the exact position of the detached shell 
with respect to Y CVn. 

We find that the surrounding medium which slows down 
the outflow is almost at rest with respect to Y CVn. Al- 
though we cannot exclude that this external medium is local 
ISM, we suggest that it might rather be made of material 
left over from an older episode of mass loss, perhaps when 
the star was on the RGB. 

Finally, we caution that the mere presence of a detached 
shell around an AGB star is not a proof that the central star 
has undergone an enhanced episode of mass loss. We stress 
the importance of acquiring high spectral resolution data 
to complement the dust continuum images obtained in the 
far-infrared. 
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